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20  ABSTRACT  (Contlnuo  on  rovotoo  oldo  II  nocoooary  and  Idontlfy  by  block  numbor) 

— Coos  Bay  Inlet,  located  on  the  south  central  coast  of  Oregon,  provides 
tidal  flow  to  two  estuary  systems:  Coos  River  to  the  north  and  South  Slough 
to  the  south.  The  upper  reaches  of  South  Slough  constitute  a National  Marine 
Sanctuary,  administered  by  the  National  Oceanographic  and  Atmospheric  Adminis- 
tration. Charleston  Harbor  is  located  at  the  entrance  into  South  Slough  and 
is  affected  by  continual  shoaling  problems  within  the  harbor  entrance  channel. 

(Continued)] 
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ABSTRACT  (Continued) 


A two-dimensional  numerical  tidal  model  (WIFM)  was  used  to  investigate  the 
tidal  hydrodynamics  of  the  inlet  complex.  Due  to  the  nature  of  the  complex 
geometry  of  the  Coos  Bay  area,  it  was  necessary  to  develop  a capability  of 
computing  the  tidal  regime  on  a variable  grid  system.  The  primary  objective 
of  this  study  was  to  apply  WIFM  to  the  Coos  Bay-South  Slough  complex  to  pre- 
dict quantitatively  the  hydrodynamics  (exclusive  of  sediment  transport  and 
wave  action)  of  the  tidal  flow  in  the  system  and  hence  draw  a comparison 
between  existing  conditions  and  alternate  improvement  plan  conditions.  In 
general,  the  improvement  plans  provide  for: 

a.  Construction  of  specified  navigation  channels. 

b.  Alternate  breakwater  extensions  west  of  the  entrance  channe 


rr 
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c.  Alternate  detached  groins  east  of  the  entrance  channel. 

Five  alternate  plans  (plan  1 and  plans  A-D)  were  tested,  and  the  results  show 
that  the  introduction  of  any  of  the  improvement  plans  as  proposed  would  not 
produce  any  detrimental  impact  to  tidal  circulation  in  South  Slough.  The 
greatest  effect  on  the  hydrodynamics  of  the  total  system  would  occur  with 
the  installation  of  plan  1.  Results  for  plans  A-D  indicate  that  t lie  combina- 
tion of  a Charleston  breakwater  extension  and  one  or  more  groins  on  the  east 
side  of  Charleston  Channel  would  best  meet  the  District's  needs  in  that  they 
permit  more  control  over  the  alignment  of  the  Charleston  Channel.  Time  his- 
tories of  tide  elevations  and  velocities  are  presented  for  selected  stations 
throughout  the  Coos  Bay-South  Slough  complex.  Sample  circulation  patterns  for 
four  instances  in  the  tidal  cycle  are  presented  for  the  verification  conditions 
and  for  each  plan.  Total  discharge  through  several  key  ranges  in  the  system  is 
graphed  as  a function  of  time.  A wave  refraction  analysis,  using  a linear  wave 
refraction  model,  was  performed  to  determine  wave  energy  levels  penetrating  the 
entrance  to  Coos  Bay  Inlet.  Results  are  tabulated  for  six  initial  deepwater 
directions  and  six  periods.  An  attempt  was  made  to  compute  wave  orthogonals 
extending  into  Coos  Bay  Inlet.  Only  those  shallow-water  waves  whose  azimuth 
is  aligned  with  the  entrance  channel  would  propagate  within  the  inlet  system. 
Results  obtained  by  the  linear  wave  refraction  code  are  very  limited  for  sites 
such  as  Coos  Bay  due  to  the  complex  channel  formation  extending  outward  from 
the  inlet  jetties  and  shoal  areas  adjacent  to  the  main  channel  as  it  progresses 
eastward  toward  Barvicw. 
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The  model  investigation  described  herein  was  authorized  by  the 
U.  S.  Army  Engineer  District,  Portland,  and  conducted  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  in  the  Wave  Dynamics  Divi- 
sion (WDD) , Hydraulics  Laboratory,  under  the  direction  of  Mr.  H.  B. 
Simmons,  Chief  of  the  Hydraulics  Laboratory,  and  Dr.  R.  W.  Whalin, 

Chief  of  the  Wave  Dynamics  Division. 

The  investigation  was  performed  and  this  report  prepared  by 
Mr.  H.  Lee  Butler,  WDD.  Mr.  George  Fisackerly,  Estuaries  Division, 
supervised  the  field  survey.  The  numerical  computations  associated 
with  this  work  were  performed  on  a Control  Data  Cyber  176  located  at 
the  Air  Force  Weapons  Laboratory  (AFWL) , Kirtland  Air  Force  Base, 

New  Mexico. 

Commander  and  Director  of  WES  during  the  course  of  the  investi- 
gation and  the  preparation  and  publication  of  this  report  was 
COL  John  L.  Cannon.  Technical  Director  was  Mr.  Fred  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply By  To  Obtain 


feet 

0.3048 

metres 

feet  per  second 

0.3048 

metres 

per  second 

square  miles  (U. 

S. 

statute) 

2.589988 

square 

kilometres 

square  feet  per 

second 

0.09290304 

square 

metres  per 

second 

feet  per  second 

per 

second 

0.3048 

metres 

per  second 

per  second 

NUMERICAL  SIMULATION  OI;  Till;  COOS  BAY-SOUTH  SLOUGH  COMPLEX 


Background 

1.  This  study,  funded  by  the  U.  S.  Army  engineer  District, 
Portland,  Engineering  Division  (NPPENJ,  was  initiated  for  the  purpose 
of  investigating  the  effect  of  several  alternative  training  dikes  and/or 
breakwaters  on  tidal  circulation  in  the  Coos  Bay-South  Slough  area. 

Coos  Bay  Inlet  is  a natural  inlet  located  on  the  south  central  coast  of 
Oregon  (Figure  1).  The  inlet  traffics  heavy  shipping  up  the  Coos  Bay 
River  to  North  Bend  and  small  craft  (commercial  fishing  and  pleasure 
craft)  to  Charleston  Harbor  just  south  of  the  inlet  entrance.  Depths 
in  the  entrance  channel  and  the  Coos  River  ranges  reach  60  ft*  NGVD 
(National  Geodetic  Vertical  Datum).  This  reference  to  depth  means  60  ft 
below  the  plane  of  NGVD  where  NGVD  is  a fixed  reference  adopted  as  a 
standard  geodetic  datum  for  heights.  The  geodetic  datum  is  fixed  and 
does  not  take  into  account  the  changing  stands  of  sea  level.  Since  the 
geodetic  datum  represents  a best  fit  over  a broad  area,  the  relation 
between  the  geodetic  datum  and  local  mean  sea  level  is  not  consistent 
from  one  location  to  another  in  either  time  or  space.  The  last  general 
adjustment  taken  into  account  for  the  data  presented  in  this  report  was 
the  1947  adjustment.  South  Slough  is  a shallow  body  of  water  connected 
to  the  system  via  a single  channel  between  the  cities  of  Charleston  on 
the  west  and  Barview  on  the  east.  The  upper  reaches  of  South  Slough 
constitute  a National  Marine  Sanctuary,  administered  by  the  National 
Oceanographic  and  Atmospheric  Administration,  and  the  first  created 
under  the  Coastal  Zone  Management  Act.  South  Slough  is  characterized 
by  narrow  channels  winding  their  way  through  mud  flats  and  sandbars. 

The  region  to  be  modeled  encompasses  an  area  of  over  46  square  miles. 

* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  3. 
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2.  Traffic  into  and  out  of  Charleston  Harbor  has  been  affected 
by  the  shifting  of  Charleston  Channel  (that  part  of  the  South  Slough 
channel  from  the  main  Coos  Bay  channel  to  Charleston  Bridge)  out  of 
project  alignment.  Continual  shoaling  problems  exist  within  the  chan- 
nel. Cost  of  additional  maintenance  may  warrant  installation  of  struc- 
tures that  would  alleviate  the  problem.  A secondary  benefit  from  such 
structures  would  be  a potential  reduction  in  wave  damage  in  the 
Charleston  boat  basin.  Such  plans  must  be  analyzed  regarding  their 
impact  on  the  tidal  hydraulics  of  the  entire  system  and  of  South  Slough 
in  particular. 

3.  A two-dimensional  numerical  tidal  model  (WIFM)  developed  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  was  used  to  per- 
form the  analyses.  Hue  to  the  nature  of  the  complex  geometry  of  the 
Coos  Bay  area,  it  was  necessary  to  develop  a capability  of  computing 
tidal  hydrodynamics  on  a variable  grid  system.  The  technique  for  this 
scheme  was  taken  from  Wanstrath  and  introduced  into  WIFM. 

Object ives 

4.  The  primary  objective  of  this  study  was  to  apply  WIFM  to  the 
Coos  Bay-South  Slougn  complex  to  predict  quantitatively  the  hydro- 
dynamics (exclusive  of  sediment  transport  and  wave  action)  of  the  tidal 
flow  in  the  system,  and  hence  draw  a comparison  between  existing  condi- 
tions and  alternate  plan  conditions.  Specifically,  the  following  tasks 
were  accomplished: 

a.  A field  survey  for  obtaining  prototype  tidal  and  velocity 
data 

b.  A numerical  investigation  of  tidal  circulation  including 
verification  and  testing  of  five  alternate  improvement 
plans 

£.  A wave  refraction  study  for  the  area  around  the  entrance 
channel 
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PART  II:  COMPUTATIONAL  TECHNIQUES 


Equations  of  Motion 


S.  The  hydrodynamic  equations  used  in  this  work  are  the  classical 
shallow-water  wave  equations.  Expressed  in  a Cartesian  coordinate  sys- 
tem (Figure  2),  the  equations  of  motion  are  given  as 


MOMENTUM 


CONTINUITY 


WATER  SURFACE 


BENCHMARK  DATUM 


Figure  2.  Coordinate  system  for  problem  formulation 


' 


The  fluid  has  been  assumed  incompressible  and  homogeneous  and  all  den- 
sity variations  are  neglected.  In  these  equations  U and  V are  the 
vertically  integrated  transports  per  unit  of  width  at  time  t in  the 
x and  y directions,  respectively;  n is  the  water-surface  elevation 
with  respect  to  the  given  datum;  d = n - h is  the  total  water  depth  at 
(x.y.t);  f is  the  Coriolis  parameter;  C is  the  Chezy  coefficient; 

Fx  and  F are  terms  representing  external  forcing  functions  such  as 
wind  effects;  g is  the  acceleration  due  to  gravity;  and  R is  a term 
representing  the  rate  at  which  additional  water  is  introduced  into  or 
taken  from  the  system  (for  example,  through  rainfall  and  evaporation).* 

A complete  derivation  of  these  equations  can  be  found  in  Wang  and 
Connor . ^ 


Numerical  Approach 

t>.  Since  obtaining  a solution  to  the  governing  nonlinear  equations 
on  a highly  complex  region  is  intractable  for  a purely  analytical  ap- 
proach, a numerical  technique  is  employed.  A finite  difference  scheme 
which  treats  the  time  dependency  implicitly  for  cost-effective  simula- 
tion is  applied,  whereby,  the  dependent  variables  of  the  centered, 
alternating-direction  procedure'5  are  the  vertically  integrated  fluid 
transports  and  surface  elevations  as  a function  of  position  and  time. 

A space-staggered  grid  scheme  is  used  and  a typical  cell  is  described  in 
figure  3.  Included  in  the  model  are  actual  bathymetry  and  topography, 

□ - FLOW/UNIT  WIDTH  IN 
X-DIRECTION  (U) 

4- FLOW/UNIT  WIDTH  IN 
Y-DIRECTION  (V) 

O-SURFACE  ELEVATIONIt,  ), 

WATER  DEPTH  (d), 

FRICTIONAL  COEFFICIENT  (C  OR  n) 


Figure  3.  Cell  definition 
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* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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variable  bottom  roughness,  inertial  forces  due  to  advcctive  and  Coriolis 

accelerations,  rainfall,  and  spatial  and  time-dependent  wind  fields. 

Flooding  of  low-lying  terrain  during  a normal  tide  cycle  is  simulated. 

The  model  is  also  capable  of  treating  subgrid  barrier  effects.  Fxposed, 

submerged,  and  overtopping  barriers  can  be  represented  within  the  grid 

system.  The  formulation  of  the  flooding  model  (entitled  WIFM--WHS  Im- 

4 

plicit  Flooding  Model)  is  presented  in  detail  by  Butler. 

Variable  Grid 


7.  The  initial  version  of  WIFM  only  permitted  the  use  of  a regular 
spaced  rectilinear  grid  system.  Due  to  the  nature  of  the  Coos  Bay-South 
Slough  complex  and  the  relatively  small  dimensions  of  plan  modifications 
to  be  tested,  it  was  necessary  to  develop  a variable  grid  procedure  to 
permit  economical  simulation  of  the  study  region.  A coordinate  trans- 
formation given  by 

x = a + botC  (4) 

where  a , b , and  c are  arbitrary  constants,  is  applied  piecewise 
for  each  axis.  One  may  consider  this  transformation  as  mapping  proto- 
type space,  discretized  with  a smoothly  varying  grid,  into  computational 
space  (a-space)  employing  a regular-spaced  grid.  The  transformation  is 
such  that  in  a-space  all  derivatives  are  centered.  By  applying  a 
smoothly  varying  grid  whose  functional  as  well  as  first  derivatives  are 
continuous,  stability  problems  usually  associated  with  variable  grid 
schemes  are  eliminated.  This  type  of  transformation  permits  simulation 
of  a complex  landscape  by  locally  increasing  grid  resolution  and/or 
aligning  coordinates  along  physical  boundaries. 

8.  The  equations  of  motion  in  a-space  can  be  written  as 


MOMFjMTIJM 


<) 


r=-  - - • 1 


3V  + ] 3 / UV  \ 

3t  Uj  3a^  \d  / 

u2  V 

• - -fr  (u2  * S)m  • 

C d 2 

(6) 

CONTINUITY 

in  + 
3t 

1 3U  1 3V 

, + _ = R 

M1  9al  U2  9a2 

(7) 

where 

P1 

= li-»  b c a(Cl-D 

3a  Til 

(8) 

p2 

= iL  = b c a(c2-^ 

3a  °2  2 2 

(9) 

Minor  modifications  of  the  finite  difference  analogs  to  Equations  1-3 
are  required  to  introduce  the  factors  and  u,  appearing  in  Equa- 
tions 5-7.  The  quantities  and  u,  define  the  stretching  of  the 

regular-spaced  computational  grid  in  a-space  (spatial  steps  of  Aa^ 
and  Aci,)  to  approximate  a study  region  in  real  or  prototype  space 
(x, y space) . 

9.  A time-share  code  has  been  designed  to  calculate  the  mapping 

defined  by  Equation  4.  Each  axis  is  partitioned  into  any  number  of 
regions  for  which  the  mapping  coefficients  are  determined.  The  code 
allows  the  user  full  control  over  the  variable  spacing  along  each  axis. 
Having  mapped  both  spatial  directions  independently,  a batch  code  plots 
the  variable  grid  in  prototype  space  on  a pen  plotter  for  use  as  an 
overlay  on  bathymetric  charts  for  ease  of  digitizing  water  depths.  In 
addition,  values  for  the  quantities  and  y.,  are  given  on  punched 

cards  for  direct  use  by  program  WIFM. 

Boundary  Conditions 

10.  Various  types  of  boundary  conditions  arc  permissible  in  the 

i present  code.  These  include: 

a_.  Open  boundaries.  Water  levels  or  flow  rates  are 
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prescribed  functions  of  location  and  time  and  are  given  as 
tabular  input  to  W1FM. 

b.  Water- land  boundaries.  Such  boundaries  lie  on  cell  faces 
and  lienee  impermeable  boundaries  are  accounted  for  by 
specifying  U = 0 or  V = 0 at  the  appropriate  cell 
face.  Low-lying  terrain  may  alternately  dry  and  flood 
within  a tidal  cycle.  Inundation  is  simulated  by  making 
the  location  of  the  land-water  boundary  a function  of  the 
current  value  of  the  total  water  depth.  By  checking  water 
level  in  adjacent  cells  relative  to  ground  elevation,  a 
determination  is  made  as  to  the  possibility  of  inundation. 
If  flooding  is  possible,  the  boundary  face  is  treated  as 
open  and  computations  for  n , U , and  V are  made  for 
that  cell.  Figure  4 depicts  a graphic  representation  of 
flood  cell  treatment.  The  drying  of  cells  is  simply  the 
inverse  process.  The  code  does  allow  excess  water  to 
drain  from  "dry  cells,"  noting  that  inundation  or  drying 
occurs  when  the  adjacent  water  level  exceeds  or  is  within 
some  small  fraction  of  a foot  (0.2  ft  in  this  application). 

£.  Subgrid  barriers.  Subgrid  barriers  are  defined  along  cell 
faces  and  are  of  three  types:  exposed,  submerged,  and 
overtopping.  Exposed  barriers  are  handled  by  simply  spec- 
ifying no- flow  conditions  across  the  appropriate  cell 
face.  Submerged  barriers  are  simulated  by  controlling  the 
flow  across  cell  faces  with  the  use  of  a time-dependent 
frictional  coefficient  (Chezy  coefficient).  Overtopping 
barrier  is  a terminology  used  to  distinguish  barriers 
which  can  be  submerged  during  one  portion  of  the  tidal 
cycle  and  totally  exposed  in  another.  Actual  overtopping 
is  controlled  by  using  a broad-crested  weir  formula^  to 
specify  the  proper  flow  rate  across  the  barrier.  Once 
the  barrier  is  submerged  (or  conversely  exposed) , pro- 
cedures as  described  above  are  followed.  Figure  5 gives 
a graphic  description  of  these  barrier  conditions. 


11 


Figure  4.  Flood  cell  treatment 
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PART  III:  TIDAL  CIRCULATION  LOR  EXISTING  CONDITIONS 


Field  Survey 

11.  A prototype  field  survey  was  performed  by  the  staff  at  WES 
during  the  period  14-22  October  1976.  Eleven  tide  gages  were  installed, 
eight  of  which  were  in  the  model  area  and  were  used  in  verifying  WIFM  to 
existing  conditions.  An  additional  gage  at  Charleston  Harbor,  monitored 
by  National  Ocean  Survey  (NOS),  was  also  used  in  the  verification  pro- 
cess. The  gages  outside  the  model  limits  were  installed  to  aid  NOS  in 
establishing  Charleston  as  a NOS  reference  station.  In  turn,  NOS  pro- 
vided technical  assistance  in  establishing  the  proper  benchmarks  for  all 
gages  installed.  A 25-hour  velocity  survey  was  carried  out  on  21-22 
October.  Current  and  salinity  measurements  were  taken  at  10  locations 
within  5 ranges.  Figure  6 depicts  all  model  tide  gage  locations  and 
velocity  range  and  station  locations.  The  velocity  survey  was  conducted 
under  a calm  sea  condition.  Since  the  calm  survey  period  is  the  only 
scenario  to  be  modeled  and  proposed  modifications  will  not  affect  wind- 
driven  circulation  in  South  Slough,  wind  effects  have  been  neglected  in 
the  study.  A concurrent  effort  by  NPPEN  to  gather  hydrographic  survey 
data  provided  definition  of  the  existing  bathymetry.  The  datum  selected 
for  WIFM  was  1947  NGVU  since  most  bathymetric  charts  provided  were  re- 
lated to  this  datum. 

12.  Tidal  gages  located  at  Cape  Arago  (gage  1)  and  Empire  (gage  0) 
were  used  to  provide  data  for  driving  the  numerical  model  at  open-water 
boundaries  (Figure  6).  Tidal  readings  at  Cape  Arago  contained  a great 
deal  of  high-frequency  components.  In  order  to  extract  meaningful  tidal 
data  from  actual  readings,  a least-squares  harmonic  (LSI!)  fit  to  eight 
days  of  available  data  was  made.  Figure  7 depicts  the  Cape  Arago  LS11 
fit  for  the  duration  of  the  velocity  survey.  The  same  process  was  also 
applied  to  the  OSU  gage  (gage  2)  and  the  fitted  data  were  used  for  com- 
parison with  model  results. 
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Figure  6.  Extent  of  model  area,  gage,  and  range  locations 
for  Coos  Bay  Inlet/South  Slough 


Figure  7.  Least-squares  harmonic  fit  to  Cape  Arago  tide  data 


Compu tational  Grid 

P 

13.  The  first  step  in  applying  WIFM  to  the  study  region  was  to 
select  a computational  grid.  A variable  mesh  for  each  direction  was 
selected  with  the  prototype  grid  spacing  ranging  from  150  ft  to  900  ft. 
The  finer  mesh  was  focused  around  the  inlet  entrance  and  Charleston 
Harbor  where  proposed  structures  to  be  tested  would  likely  be  placed. 
Figure  8 shows  the  computing  grid  overlaid  on  an  area  map.  The  grid 
consists  of  13,130  incongruent  cells  having  dimensions  of  130  x 101. 
Recall  that  a coordinate  transformation  in  the  form  of  a piecewise 
"exponential  squeeze"  is  applied  independently  to  each  axial  direction 
to  map  the  variable  grid  into  an  equally  spaced  computing  grid  used 
directly  by  WIFM. 

14.  Fourteen  bathymetry  maps  of  the  model  area  were  reduced  to  the 
same  scale  (l:b000)  and  pieced  together  to  form  three  maps:  Empire  to 
Charleston  Harbor,  Charleston  Harbor  to  the  south  end  of  South  Slough, 
and  the  open  ocean  area.  Grid  overlays  for  digitizing  water  depths  were 
generated  and  plotted  on  a Calcomp  drum  plotter.  Depth  and  keyed  bottom 
characteristics  (for  defining  associated  Manning's  n)  were  tabulated  for 
each  grid  cell. 

15.  As  mentioned  previously,  the  open  ocean  boundary  was  forced 
with  the  tide  history  obtained  at  Cape  Arago,  and  the  open  water  bound- 
ary near  Empire  was  forced  with  the  tide  history  taken  at  a nearby  dock. 
All  other  boundaries  were  taken  as  movable  land-water  boundaries.  The 
effects  of  small  channels  entering  the  southernmost  end  of  South  Slough 
were  neglected  since  their  effect  on  comparing  alternate  plan  conditions 
with  existing  conditions  would  be  small. 

lb.  Various  subgrid  anomalies  exist  throughout  the  system,  such  as 
boat  docks,  loading  docks,  narrow  shoals,  submerged  jetties,  exposed 
jetties,  and  bridges.  These  features  arc  treated  as  barriers  described 
by  different  frictional  coefficients  to  simulate  their  effect  on 
circulation. 
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figure  S.  Computational  grid  for  the  Coos  Ray  model 
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Model  Verification 


17.  The  objective  of  the  verification  phase  of  any  model  study  is 
to  demonstrate  the  model's  ability  to  produce  results  that  agree  with 
known  data  for  existing  conditions.  Conditions  for  the  model  area  from 
0800  PST  on  21  October  to  0800  PST  on  22  October  were  simulated  to  es- 
tablish the  model's  predictive  capability  as  well  as  base  conditions  for 
comparing  effects  of  proposed  improvement  plans.  Comparisons  between 
prototype  and  model  were  made  for  tidal  elevations  at  9 locations  and 
current  velocities  at  10  stations.  Prototype  current  measurements  were 
taken  at  three  depths:  surface,  middepth,  and  bottom  as  exemplified  in 
Figure  9.  These  field  data  all  indicate  that  neglecting  velocity  varia- 
tions in  the  vertical  direction  (an  assumption  used  in  formulation  of 
the  numerical  model)  is  reasonable  and  justified  for  the  study  area. 

For  ease  of  presentation,  comparison  between  model  and  prototype  veloc- 
ities will  be  made  by  using  the  prototype  middepth  velocity. 

18.  Surface  elevation  results  are  depicted  in  Plates  1-3  for  gages 
0-8.  Prototype  data  for  both  tide  and  velocity  comparisons  are  depicted 
by  a solid  line  and  model  results  by  a line  annotated  with  the  symbol  x. 
All  data  were  plotted  at  half-hour  intervals.  As  mentioned  previously, 
gages  0 and  1 were  used  to  drive  the  numerical  model  and  therefore  show 
excellent  agreement.  Note  that  the  Cape  Arago  gage  (gage  1)  is  well  re- 
moved from  the  model  boundary.  Results  indicate  that  its  use  as  a driv- 
ing gage  is  a reasonable  assumption.  All  tide  gages  outside  of  South 
Slough  show  very  good  agreement  between  model  and  prototype.  Gages 
within  South  Slough  (5-8)  depict  a good  comparison  although  an  overshoot 
in  the  flood  phase  is  noted.  This  is  probably  due  to  the  cruder  grid 
representation  used  in  the  slough  and  the  inability  to  define  mud  flats 
and  other  storage  areas  accurately  due  to  lack  of  prototype  data.  Gage 
8 is  representative  of  the  hydrograph  at  llinch  Bridge.  Because  of  the 
large  size  of  computational  cells  at  the  south  end  of  the  slough,  gage  8 
cannot  be  simulated  as  accurately  as  other  gages.  The  variable  computa- 
tion grid  was  designed  to  give  finer  resolution  in  the  area  of  proposed 
structural  changes.  Fine  resolution  throughout  the  study  area  would 
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RANGE  2 STATION  E 


Figure  9.  Example  of  the  vertical  variation  in  prototype  current  measurements 


have  resulted  in  a much  larger  number  of  grid  points  forcing  WIFM  stor- 
age requirements  to  exceed  the  capacity  of  the  Cyber  176,  the  fastest 
and  largest  computer  currently  available  to  the  Corps.  The  grid  se- 
lected permits  an  economical  and  accurate  assessment  of  plan  effects 
relative  to  base  conditions. 

19.  Velocity  comparisons  for  ranges  1-5  are  presented  in  Plates  4- 
7.  It  should  be  noted  that  the  numerical  model  began  computing  from  a 
quiescent  state  at  0600  PST.  It  takes  about  four  hours  for  the  model 
to  spin  up.  This  fact  would  explain  the  poorer  comparison  in  the  first 
hours  of  the  velocity  survey.  Some  of  the  prototype  data  are  missing 
due  to  operational  problems  in  rough  seas.  Comparison  with  prototype 
velocity  data  for  most  ranges  and  stations  is  good  (recall  that  the 
middepth  reading  at  all  stations  is  used  for  a comparison  basis). 

Range  1 extends  from  Coos  Head  to  the  north  jetty  across  the  inlet 
throat.  Range  2 is  located  between  Pigeon  Point  and  North  Spit  across 
Coos  Bay.  The  low  flood  velocity  in  the  second  flood  phase  at  sta  2D 
is  an  anomaly  which  may  be  due  to  small  eddies  not  simulated  by  WIFM. 
Range  3 is  located  east  of  the  entrance  channel  to  Charleston  Harbor. 
Again,  this  area  is  prone  to  small  eddies  as  seen  in  the  prototype  data 
at  the  beginning  of  each  ebb  phase.  Range  4 is  located  near  the  en- 
trance to  South  Slough.  Range  5 is  located  near  Copes  Dock  in  South 
Slough.  For  Range  5 the  numerical  model  results  show  a weak  flood  vel- 
ocity and  a stronger  ebb  velocity.  Sta  5J  is  located  at  the  entrance  to 
a small  lagoon  in  front  of  Copes  Dock  and  sta  5K  is  located  in  the  main 
slough  channel  200  ft  west  of  sta  5J.  There  is  no  indication  that  a 
freshwater  inflow  from  the  south  or  an  anomalous  wind  condition  caused 
the  large  ebb  velocity.  The  difference  in  tidal  range  between  the  first 
and  second  cycle  appears  insufficient  to  cause  the  prototype  ebb  veloc- 
ity in  the  first  phase  whose  peak  magnitude  is  four  times  that  experi- 
enced in  the  second  ebb  phase.  Although  some  differences  between  proto- 
type data  and  numerical  model  predictions  occurred,  the  overall  quality 
of  the  model  verification  was  good.  The  greater  differences  occurred  in 
the  areas  of  cruder  spatial  resolution.  As  long  as  changes  in  hydrau- 
lics caused  by  the  introduction  of  control  structures  do  not  extend  into 
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these  areas,  the  verification  results  can  be  used,  with  a high  degree  of 
reliability,  as  base  conditions  for  plan  effect  comparisons . Hence,  for 
the  remainder  of  the  report,  verification  results  will  be  referred  to  as 
base  conditions. 

20.  Additional  velocity  stations  were  selected  for  comparing  plan 
with  base  conditions.  Results  will  be  presented  in  PART  IV  of  this 
report.  Various  ranges  for  comparison  of  volumetric  discharges  also 
were  plotted  for  use  in  analyzing  plan  effects. 

21.  l-'or  plotting  circulation  patterns,  the  computational  region 
was  subdivided  into  three  areas:  Coos  Bay  River  Reach  (extending  from 
Lmpire  to  Charleston  Harbor),  South  Slough  (from  Charleston  Harbor  to 
flinch  Bridge),  and  an  enlargement  of  the  Inlet/Charleston  Harbor  area. 
Two  types  of  flow  field  plots  were  made:  the  first  displaying  transport 
vectors  and  the  second  displaying  velocity  vectors  (similar  to  surface 
current  patterns).  Plots  for  each  simulation  hour  were  produced  on 
microfiche,  55  mm  color  film,  and  black  and  white  film.  Plates  8-10 
display  example  circulation  patterns  (transport  and  velocity)  for  the 
inlet  subregion  at  three  time  instances  in  the  tidal  cycle:  ebb,  flood, 
and  near  slack  water.  These  plots  indicate  the  amount  and  direction  of 
water  flowing  through  a given  cell  cross  section  at  a given  time.  A 
vector  with  length  equal  to  a minimum  cell  width  represents  a flow  per 
unit  width  of  100  ft'/sec  (transport  plots)  and  a velocity  of  3 ft/sec 
(velocity  plots).  A complete  set  of  circulation  plots  is  on  file  at 
both  WHS  and  NPPEN.  A few  geographic  landmarks  have  been  identified  in 
Plate  8 to  aid  in  plot  orientation.  Heavy  straight  lines  on  each  cir- 
culation plot  depict  existing  or  proposed  structures  (jetties,  break- 
waters, bridges,  boat  docks,  etc.).  Such  features  shown  in  Plates  8-10 
are  all  existing  structures. 


PART  IV: 


ALTERNATE  IMPROVEMENT  PLANS 


22.  For  the  purpose  of  reducing  maintenance  costs  associated  with 
the  navigation  channel  into  Charleston  Harbor,  NPPEN  has  proposed  test- 
ing five  alternate  improvement  plans.  Each  plan  includes  dredging  a new 
straight  channel  to  the  harbor  as  well  as  some  combination  of  break- 
waters and/or  groins.  Figure  10  displays  all  plan  conditions  and  major 
velocity  stations  used  in  comparing  plan  with  base  conditions.  The  five 
alternative  plans  are  summarized  as  follows: 

a.  1’  1 an  1 --  An  L-shaped  breakwater  extending  from  land  at  a 
point  approximately  1,000  ft  northward  from  the  westward 
end  of  the  existing  breakwater  toward  Fossil  Point  for 
2,600  ft  and  then  north  at  a right  angle  toward  North  Spit 
for  1,000  ft  and  three  east-west  groins  (approximately  400 
ft  long)  located  between  the  existing  submerged  jetty  and 
the  entrance  to  Charleston  Harbor  on  the  east  side  of 
Charleston  Channel. 

1}.  Plan  A --  A 1,000-ft  extension  of  the  existing  Charleston 
breakwater  toward  the  north. 

£.  Plan  B --  The  breakwater  extension  in  plan  A and  a 400-ft 
east-west  groin  located  on  the  east  side  of  the  entrance 
channel  at  the  end  of  the  new  extension. 

d.  Plan  C --  Placement  of  a groin  as  described  in  plan  B 
without  the  breakwater  extension. 

£.  Plan  D --  Conditions  as  for  plan  B with  the  addition 
of  another  400-ft  east-west  groin  at  the  entrance  of 
Charleston  Harbor. 

The  new  channel  for  plan  1 would  be  dredged  to  -15  ft  NGVD,  and  to 
-14  ft  NGVD  for  all  other  plans. 

23.  Plates  11-19  show  the  comparison  of  surface  elevations  with 
plans  installed  using  verification  results  as  a base  condition.  No 
significant  changes  are  noted  between  plan  and  base  conditions. 

24.  Plates  20-29  depict  the  comparison  between  plan  and  base  con- 
ditions for  current  velocities  at  base  station  locations.  Again,  no 
significant  changes  are  noted  with  any  plan  installed,  with  the  excep- 
tion of  sta  3G.  This  station  is  located  between  or  near  groins  in 
plans  1,  B,  C,  and  0.  As  expected,  the  velocity  is  reduced  with  the 
installation  of  groins  east  of  the  new  channel. 
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Figure  10.  Proposed  plans  for  Coos  Bay  Inlet  with 
major  velocity  stations  depicted 


J5.  Additional  velocity  stations  were  selected  for  comparing  plan 

with  base  conditions.  Station  locations  (I* IV)  are  shown  in 

figure  b and  results  are  depicted  in  Plates  30-37.  Sta  P is  located  at 
the  tip  of  the  1,-shaped  breakwater  in  plan  1 and  the  velocity  at  sta  P 
is  reduced  by  introducing  the  plan.  Since  the  breakwater  extension  in 
plans  A,  B,  and  0 does  not  reach  the  area  of  sta  P,  little  change  is 
noted  with  t he  introduction  of  those  plans.  Sta  Q and  R arc  located 
between  or  near  training  dikes  in  all  plans  except  plan  A.  Results  in- 
dicate t he  suppression  of  flow  in  these  areas.  Sta  S is  located  at  t he 
confluence  of  the  main  channel  and  the  Charleston  Harbor  (Cll)  channel. 
Little  change  is  noted  with  the  introduction  of  plans  AD.  However, 
plan  1 greatly  increases  the  local  velocity  and  discharge  into  and  out 
of  the  Cll  channel  since  the  tip  of  the  I, -shaped  breakwater  nearly 
reaches  this  junction.  Sta  T is  located  in  the  till  channel  adjacent  to 
the  end  of  t lie  present  Charleston  breakwater.  Introduction  of  a groin 
east  of  the  channel  will  greatly  increase  the  current  velocity  in  the 
channel  as  noted  in  plans  1 and  l).  Sta  II  is  at  the  entrance  to  South 
Slough  and  sta  V at  the  entrance  to  doe  Ney  Slough.  No  significant 
velocity  change  is  noted  with  the  introduction  of  any  of  the  proposed 
plans.  Sta  W is  located  in  the  Cll  channel  at  the  entrance  to  the  boat 
basin.  Introduction  of  a groin  opposite  the  end  of  the  present 
Charleston  breakwater  slightly  increases  current  velocity  into  and  out 
of  the  boat  basin. 

Plates  38-(>7  display  vector  flow  patterns  (transport  and  vel 
oeity)  at  three  time  instances  in  the  tidal  cycle  for  each  plan.  Pro 
posed  structures  are  depicted  by  additional  heavy  straight  lines  in 
comparison  with  Plates  S 10.  Also  given  are  difference  plots  (traits 
port  and  velocity)  for  each  plan.  These  differences  were  formed  bv 
subtracting  vector  components  of  base  calculations  from  plan  results. 
The  direction  of  the  resulting  vector  will  be  reversed  relative  to  the 
proper  flow  direction  if  the  flow  magnitude  for  the  base  condition  is 
greater  than  that  for  plan  conditions.  All  difference  plots  indicate 
that  the  impact  on  circulation  due  to  plan  installment  i s confined  to 
a localised  area  near  the  structures. 
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27.  Several  discharge  ranges  were  selected  to  present  graphs  of 
volumetric  discharge  as  a function  of  time.  Plates  68-72  depict  the 
discharge  through  five  ranges  (R3-R6,  R9)  whose  locations  are  given  in 
Figure  6.  No  change  is  noted  with  the  installation  of  any  plan  for 
R3  and  R9  (inlet  throat  and  Coos  Bay-Empire  ranges,  respectively). 

Ranges  R4,  R5,  and  R6  are  located  at  the  entrance  to  Charleston  Harbor, 
South  Slough,  and  Joe  Ney  Slough,  respectively.  Little  change  in  total 
discharge  is  noted  for  these  ranges,  with  plan  1 showing  the  largest  im- 
pact at  any  one  instant  in  time.  Results  from  these  and  other  ranges 
indicate  that  the  tidal  prism  in  South  Slough  will  not  be  altered  by 
construction  of  any  of  the  plans  tested. 
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PART  V:  WAVE  REFRACTION  ANALYSIS 


Methodology 

28.  A wave  refraction  analysis  was  performed  to  determine  wave 
energy  levels  penetrating  the  entrance  to  Coos  Bay  Inlet.  This  phase  of 
the  study  was  conducted  using  a linear  wave  refraction  program  developed 
by  Dobson. Effects  of  both  reflection  and  diffraction  are  neglected, 
and  thus  it  is  assumed  that  no  energy  is  transmitted  along  the  wave 
crest.  These  assumptions  lead  to  the  wave  height  being  defined  by  the 
relat ion 


H = H K K 
o s r 

where  Hq  is  the  wave  height  in  deep  water,  Kg  is  the  shoaling  co- 
efficient, and  Kr  is  the  refraction  coefficient.  Refraction  diagrams 
were  produced  from  two  grids  as  depicted  in  Plate  73.  The  open  ocean 
grid  extends  seaward  to  water  depths  of  600  ft  and  the  inlet  entrance 
grid  extended  from  the  seaward  end  of  the  jetties  forming  the  inlet 
entrance  to  Barview.  The  inlet  grid  was  developed  to  compute  the  dis- 
position of  wave  energy  penetrating  the  inlet. 

Deepwater  Grid  Results 

29.  Wave  orthogonals  were  calculated  and  plotted  from  initial 
deepwater  directions  SW,  WSW,  W,  WNW,  NW,  and  NNW,  and  for  periods  of 
8,  10,  12,  14,  16,  and  18  sec.  Figures  for  the  deepwater  refraction 
patterns  were  assembled  in  a separate  report  to  NPPEN.  Plates  74-79 
display  samples  from  each  deepwater  direction.  Table  1 summarizes  the 
results  for  this  phase  of  the  analysis.  Values  given  in  the  table  are 
representative  of  an  average  of  orthogonals  in  the  harbor  entrance. 

The  wave-height  adjustment  factor  is  the  product  of  and  , and 

can  be  applied  to  any  deepwater  wave  height  to  obtain  the  corresponding 
shallow-water  wave  height.  Since  linear  theory  breaks  down  in 
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convergence  areas,  a maximum  refraction  coefficient  of  1.4  and  a minimum 
of  0.4  should  be  taken  as  reasonable  limiting  values.'7 


- 


Inlet  Grid  Results 

30.  For  four  deepwater  directions,  SW,  WSW,  NW,  and  NNW,  the  di- 
rection of  the  shallow-water  wave  at  the  inlet  opening  was  at  such  an 
angle  that  the  computer  code  could  not  compute  orthogonals  extending 
much  beyond  Coos  Head  just  inside  the  inlet  entrance.  Four  examples  are 
depicted  in  Plates  80-83  for  a period  of  12  sec.  The  wave  front  azimuth 
at  the  inlet  opening  was  computed  by  averaging  the  shallow-water  azi- 
muths for  each  deepwater  direction.  The  shallow-water  azimuths  used  in 
the  inlet  grid  runs  are  tabulated  below. 

Shallow-Water 


Deepwater  Direction 

Azimuth 

SW  (225°) 

272.5° 

WSW  (247.5°) 

276.5° 

Average  of  W and  WNW 
(270°  and  292.5°) 

297.5° 

NW  (315°) 

315.5° 

NNW  (337.5°) 

321.5° 

31.  Wave  energy  approaching  the  interior  of  the  bay  is  concen- 
trated from  directions  W to  WNW  (75  percent)  as  given  by  National  Marine 

g 

Consultants.  The  average  shallow-water  azimuth  resulting  from  these 
deepwater  directions  is  297.5°.  Since  shallow-water  waves  from  such  an 
azimuth  begin  propagating  directly  up  the  entrance  channel,  these  condi- 
tions were  considered  the  ones  in  which  maximum  wave  energy  would  reach 
the  Barview  area.  Refraction  diagrams  for  wave  periods  from  8-18  sec 
are  given  in  Plates  84-90.  Note  that  two  figures  are  presented  for  the 
14-sec  period.  The  second  diagram  (Plate  88)  indicates  the  sensitivity 
in  choosing  the  spacing  between  wave  orthogonals  at  the  beginning  of  the 
computation.  Since  each  wave  ray  is  computed  independently,  it  is  pos- 
sible to  compute  crossovers  as  exhibited.  One  can  consider  each  ray  as 
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representing  the  movement  of  a "packet"  of  energy,  and  thus  it  is  not 
inconsistent  for  such  crossovers  to  occur.  The  results  obtained  by  the 
linear  wave  refraction  code  are  very  limited  for  sites  such  as  Coos  Bay 
due  to  the  complex  channel  formation  extending  outward  from  the  inlet 
jetties  and  shoal  areas  adjacent  to  the  main  channel  as  it  progresses 
toward  Barview. 
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PART  VI:  CONCLUSIONS 


32.  The  Coos  Bay-South  Slough  numerical  model  was  verified  to  re- 
produce satisfactorily  the  hydrodynamic  response  of  the  Coos  Bay  inlet 
system  to  a specified  astronomical  tide.  The  verification  procedure 
substantiates  the  model's  ability  to  reliably  predict  the  impact  of 
proposed  modifications  to  the  Charleston  Channel  area  on  tidal  eleva- 
tions, current  velocities,  and  circulation  within  the  study  region. 

33.  Results  of  the  five  plans  tested  show  that  construction  of 
any  of  these  improvement  plans  tested  will  have  a minimal  effect  on 
tidal  circulation  in  South  Slough;  most  velocity  and  circulation  changes 
concerned  with  these  plans  were  confined  to  the  local  vicinity  of  the 
proposed  modifications.  The  greatest  effect  on  the  hydrodynamics  of  the 
total  system  would  occur  with  construction  of  plan  1.  Results  for  plans 
A-D  indicate  that  the  combination  of  a Charleston  breakwater  extension 
and  one  or  more  groins  on  the  cast  side  of  Charleston  Channel  would  best 
meet  the  District's  needs.  These  plans  are  less  expensive  to  construct 
than  plan  1 yet  still  permit  control  over  alignment  of  the  Charleston 
Channc 1 . 

34.  No  conclusions  regarding  reduction  of  shoaling  problems  in  the 
Charleston  Channel  or  of  wave  damage  in  the  Charleston  boat  basin  can  be 
made.  The  tidal  circulation  model  is  not  designed  to  address  these 
problems  and  wave  refraction  models  are  severely  limited  by  inherent 
assumptions  of  no  energy  transfer  across  wave  orthogonals.  Consequently 
results  from  wave  refraction  analyses  are  of  limited  value  in  quanti- 
tatively evaluating  the  effect  of  various  plans  on  wave  energy  in 
Charleston  boat  basin.  However,  inspection  of  the  circulation  patterns 
and  wave  refraction  analyses  coupled  with  experienced  coastal  engineer- 
ing judgment  may  lend  some  insight  into  these  problems. 
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Table  1 

Summary  of  Refraction  and  Shoaling  Analysis  at  the 
Entrance  of  Coos  Bay  Harbor,  Oregon 


Deepwater 

Direction 

Wave 

Period 

sec 

Shallow-Water 

Azimuth 

deg 

Refraction 

Coefficient 

Shoaling 

Coefficient 

Wave-Height 

Adjustment 

Factor 

SW(225°) 

8 

268 

0.48 

0.95 

0.46 

10 

263 

0.30* 

1.01 

0.30 

12 

267 

0.26* 

1.08 

0.28 

14 

279 

0.58 

1.09 

0.63 

16 

280 

0.67 

1.11 

0.74 

18 

278 

0.66 

1.34 

0.88 

i WSW (247.5°) 

8 

266 

0.38* 

0.96 

0.37 

10 

270 

0.58 

0.96 

0.56 

12 

274 

0.73 

1.08 

0.79 

1 

14 

278 

0.66 

1.14 

0.75 

16 

282 

0.61 

1.17 

0.71 

18 

289 

0.48 

1.17 

0.56 

WEST(270°) 

8 

281 

0.86 

0.97 

0.83 

10 

288 

0.62 

0.98 

0.61 

12 

291 

0.68 

1.02 

0.69 

14 

295 

0.46 

1.10 

0.51 

16 

297 

0.41 

1.11 

0.46 

18 

299 

0.80 

1.22 

0.98 

WNW (292 . 5°) 

8 

302 

0.71 

0.93 

0.66 

10 

303 

0.62 

0.97 

0.60 

12 

300 

1.01 

1.03 

1.04 

14 

301 

0.77 

1.08 

0.83 

16 

305 

0.92 

1.11 

1.02 

18 

307 

1.10 

1.16 

1.28 

NW(31S°) 

8 

315 

0.84 

0.95 

0.80 

10 

313 

0.70 

1.00 

0.70 

12 

315 

0.76 

1.03 

0.72 

14 

316 

0.58 

1.05 

0.61 

16 

317 

0.72 

1.18 

0.85 

18 

317 

0.52 

1.19 

0.62 

NNW (337. 5°) 

8 

323 

0.97 

0.94 

0.91 

10 

327 

0.81 

1.00 

0.81 

12 

322 

0.71 

1.01 

0.72 

14 

318 

0.62 

1.10 

0.68 

16 

320 

0.56 

1.13 

0.63 

18 

319 

0.61 

1.20 

0.73 

* A minimum  refraction  coefficient  of  0.40  should  be  used. 
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APPENDIX  A:  NOTATION 

a,b,c 

Arbitrary  constants 

b 

Barrier  height  (Figure  5) 

C 

Chezy  frictional  coefficient 

1 C° 

Admittance  coefficient 

d 

Total  water  depth 

dll 

Water  depth  over  barrier 

f 

Coriolis  parameter 

F F 
rx»ry 

External  forcing  functions 

F F 

V “2 

External  forcing  functions  in  computational  space 

g 

Acceleration  due  to  gravity 

h 

Land-surface  elevation 

Ho 

Wave  height  in  deep  water 

Kr 

Refraction  coefficient 

Ks 

Shoaling  coefficient 

M,N 

Indices  denoting  spatial  increments  in  the  x and  y direction 

n 

Manning's  frictional  coefficient 

R 

Rate  at  which  additional  water  is  introduced  into  or  taken 
from  the  system  (rainfall  minus  evaporation) 

t 

Time 

u,v 

Vertically  integrated  horizontal  transports  per  unit  width 

W 

Vertical  transport 

x,y,z 

Cartesian  coordinates 

a 

Independent  variable  in  computational  space 

A , A 

V a2 

Spatial  steps  in  computational  space 

e 

Small  depth  of  water 

n 

Water-surface  elevation 

Uj,U2 

I 

Grid  transformation  stretching  factors 
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